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ABSTRACT 


The purpose of this study was to predict if subcoded cryogenic liquid entering the bottom ctf a 
storage tank will destroy the th^mal stratification of the tank. 

After an extensive literature search, a fcamula for maximum oitical Reynolds Number whidi used 
to predict the destratification of a cryogenic tank was found. Example calculations and graphics 

to detomine the mixing of fluid in the tank were presented. 
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SUMMARY 

1 . Formula erf" Critical Reynolds number to predict destratificaticgi was shown. 

2. Gratis of Reynolds number inlet flow pipe, and critical Reynolds number versus 
different inlet pipe diameter were presented. 

3. A sample of calculation of liquid hydrogen tank, and conclusion was given. 
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1. INTRODUCTION 


1.1 BACKGROUND 

Thermal stratification exists in cryogenic fluid in a storage tank when there is a temperature 
differential between fluids on the top and the bottom of the storage tank. This temperature 
differential is caused by the buoyant forces resulting from the density difference between hot and 
cold fluid. Because of the buoyant forces, it is evident that some degree of stratification exist in 
almost all ayogenic stcvage tanks unless there is a mixing force. 

This study invcrfved a cryogenic refiigouticm subcooling systems which take liquid frtxn die top of 
a cryogenic storage tank and passed it through a refrigerator and back into the bottom of the tank. 
Figure 1 shows a schematic of piping connecting the storage tank. In order to make the 
refrigeration process more efficient the system should provide as high a temperature difference 
between the cryogen at the inlet and outlet of the refrigerator as possible. Mixing of the fluid in 
the tank would lower the temperature of the cryogen going into the refrigerator and lower the 
temperature difference. This study was made to determine how to analyze the effect of subcooled 
liquid entering the bottom erf* the tank on thermal stratification in the tank. 
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Fig. 1 Schematic of Cryogenic Tank System 
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1.2 OBJECTIVE 

The object of this study was to gather applicable engineering information and data and see if 
cryogen entering the bottom of a cryogenic storage tank, will cause the hot and cold liquid to mix. 
The cryogenic liquids studied were liquid oxygen and liquid hydrogen. Typical tank sizes were 
sdected. 


2. CALCULATIONS 


2. 1 CRITICAL REYNOLDS NUMBER AND REYNOLDS NUMBER 


A critical Reynolds number of the inlet flow pipe was defined in reference [3] as that value where 
the buoyant force of tank liquid becomes unimportant For a Reynolds number greater than this 
critical value, the system may be assumed to mix completely. Below this critical Reynolds 
number, the buoyant forces may be strong enough to limit the degree of mixing. The critical 
Reynolds numb^ was given by a correlation available in reference [3]. 


The critical Reynolds number is given by 


Re* - 0.912^— r Pt ^ Gr « 




r -2 


( 1 ) 


where Re = Critical Reyndds number of inlet pipe 

Pr = Prandtl number 
D = Tank diameter 
H = Tank height 
Gr* = Modified Grashof number 


The modified Grashof number is given in reference [3] as 


Gr*~ 


kv^ 


( 2 ) 
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where: g = Gravitational acceleration 

P - Coefficient of thermal expansion ctf liquid 
q” = Heat flux 
L = Liquid hdght 
k - Thermal ccMiductivity of liquid 
V = Kinematic viscosity of liquid 


The Reynolds number is given as 


where 



(3) 


V = Flow velocity 


d = Inside diameter of inlet jripe 
V «* Kinematic viscosity 


2.2 SAMPLE OF CAIjCULATION 


2.2. 1 THE CRITICAL REYNOLDS NUMBER FOR LIQUID HYDROGEN 

The critical Reyncrfds number for liquid hydrogen is calculated from data given in Table 1. 
Pressme, temperature, and heat flux are assumed. 


TABLE 1 
Term 

Svmbol 

LH2 Tank 
Unit 

Value 

Pressure in tank 

P 

psia 

35 

Outlet temperature 

To 

R 

27 

Tank diameter 

D 

ft 

10 

Tank height 

H 

ft 

37.75 

Tank vcdume 

Vt 

ft^ 

2669 

TABLE 1 (Continued) 




Reference 
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Symbol Unit 


Term 

Inside dia. of inlet pipe. 
Inlet flow rate 
Thermal expansion 
Density 
Specific heat 
Dynamic viscosity 
Thermal conductivity 
Prandtl number 
Heat flux 

Kinematic viscosity 
Tank liquid height 
Tank liquid volume 
Gravitation acceleration 


d 

ft 

Q 

f^ /hr. 

P 

1/R 

P 

Ibm/ft? 

Cp 

Btu/lbm. R 

p 

Ibm/ft hr. 

k 

Btu/fthr.R 

Pr 


q- 

BtuAir. 

V 

ft.2/hr. 

L 

ft 

VI 

fP 

g 

ft/hr^ 


Value Reference 


0.5 

1299.5 

0.006256 

[1] 

4.758 

[1] 

1.6595 

[1] 

0.0539 

[1] 

0.04795 

[1] 

1.867 

[1] 

14.23 

0.01133 

[1] 

34.5 

2562 

417312,000 



Frcxn Equation (2) 


Gr* 


kv^ 


4.173(l0 *)^(6.256Xl0)~*(^)(34.5)'';i‘(14.23)-^ 


Gr* = 8.549 x 10^^ 


From Equation (1) 


Re* 


Re - 


0.9121 




2 

Gr*~^ 


0.912(^) ’(1.867)f[8.549(10)*]5 


7 


Re (Critical Reynolds number) = 1.0256 x 10^ 

(for the selected tank with D = 10 ft, H = 37.75 ft and L = 34.5 ft) 


577 



10 


2.2.2 THE CRITICAL REYNOLDS NUMBER AND REYNOLDS NUMBER FOR UQUED OXYGEN 

The calculation for Critical Reynolds number and Reynolds numbers at different inlet pipe 
diameters are shown in the Appendix section. 

3. RESULTS AND CONCLUSION 

3.1 RESULTS 

The Critical Reynolds number and Reynolds number for both liquid hydrogen and liquid oxygen 
are calculated by using the Microsoft Excel spreadsheet are completed (see Appendix). Graphs of 
critical Reynolds number and Reynolds number versus inlet pipe diameter for both LH2 and L02 
are shown in Hg. 2 and Hg. 3. 



INSIDE DIAMETER OF INLET PIPE (IN.) 
Fig. 2 Liquid Hydrogen Tank Mixing Curve 


578 




REYNOLDS NUMBER 



Fig. 3 Liquid Oxygen Tank Mixing Curve 


579 


12 


3.2 CONCLUSION 

In conclusicxi, form Fig. 2, it indicates that for liquid hydrogen tank an inlet pipe diameter about 
1.8 in. or less, the tank fluid will have a complete mix. Diameter of die inlet pipe greater the 1.8 in 
will be large mough to prevoat axn|:dete mixing. 

Hg. 3 indicates that for liquid oxygen tank an inlet pipe diameter of about 4 in. or less, the tank 
fluid will have a complete mix. Diameter of the inlet pipe greater the 4 in. will prevent complete 
mixing. 

The dimensions of the selected hydrogen tank are D = 10 ft, H = 37.75 ft, and L = 34.5 ft 
The pressure and temperature at the tank inlet were assumed to be 25 psia and 27 R respectively. 
The dimensions of the selected oxygen tank are D = 10 ft., H = 23.83 ft, and L = 21.22 ft 
The pressure and temperature at the tank inlet were assumed to be 35 psia and 141 R respectively. 
Results shown in this report are used as a guide to place a reasonable predicticm only. 
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CALCULATION OF MODIFIED GRASHOF NUMBER 


Formula: 


GR* 


kv^ 


CALCULATION OF CRITICAL REYNOLDS NUMBER 
Formula: 

/ rv I 2 

Re* - 0*912^^ Rr 3 Gr *3 

TABLE 1 LIQUID HYDROGEN TANK 


Term 

Symbol 

Unit 

Value 

Referen 

Pressure 

P 

psia 

35 


Outlet temperature 

To 

R 

27 


Tank diameter 

D 

ft 

10 


Tank height 

H 

ft 

37.75 


Tank volume 

V 

CufL 

2669 


Inside dia. o( inlet pipe. 

d 

ft 

0.5 


Inlet flow rate 

Q 

cu fL/hr. 

1299.5 


Thermal expansicHi 
coefficient 

P 

1/R 

0.006256 

[1] 

Doisity 

P 

Ibm/cu ft 

4.758 

[1] 

Spedflc heat 

Cp 

Btu/lbm. R 

1.6595 

[1] 

Dynamic viscosity 

p 

Ibm/ft. hr. 

0.0539 

[1] 

Thermal conductivity 

k 

Btu/fLhr.R 

0.04795 

[1] 

Prandtl number 

Pr 


1.867 

[1] 

Heat flux 

q” 

Btu/hr. sq ft 

14.23 


Kinematic viscosity 

V 

sq ft/hr. 

0.01133 

[1] 

Tank liquid height 

L 

ft 

34.5 


Tank liquid vcrfume 

Vt 

cu ft 

2562 


Giavitaticxi acceleration 

g 

ft/sqhr 

417312,000 
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LH2 Tank 


Modified Grashof number 
Gr* = 8.549 x 10 l 8 


Critical Reynolds number of the LH2 tank 
with selected D, H, and L 

Re = 1.0256 x 10<» 


Data for curve in Fig. 2 


Inside Pipe 
Diameter 
(in.) 

Area (sq ft) 

velocity (ft/s) 

Reynolds No. 

Critical Re No. 

1 

5.45E-03 

6.62E+01 

1.75E+06 

1.03E+06 

2 

2.18E-02 

1.65E+01 

8.76E+05 

1.03E+06 

3 

4.91E-02 

7.35E+00 

5.84E+05 

1.03E+06 

4 

8.73E-02 

4.14E+00 

4.38E+05 

1.Q3E+06 

5 

1.36E-01 

2.65E+00 

3.50E+05 

1.03E+06 

6 

1.96E-01 

1.84E+00 

2.92E+05 

1.03E+06 

7 

2.67E-01 

1.35E+00 

2.50E+05 

1.03E+06 

8 

3.49E-01 

1.03E+00 

2.19E+05 

1.03E+06 

9 

4.42E-01 

8.17E-01 

1.95E+05 

1.03E+06 

10 

5.45E-01 

6.62E-01 

1.75E+05 

1.03E+06 

11 

6.60E-01 

5.47E-01 

1.59E+05 

1.Q3E+06 

12 

7.85E-01 

4.60E-01 

1.46E+05 

1.03E+06 
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Term 

Pressure 

Outlet tempoature 
Tank diameter 
Tank height 
Tank volume 
Inside dia. of inlet pipe 
Tank inlet flow rate 
Thoroal expansion 

coefficient 

Daisity 
Specific heat 

Dynamic viscosity 
Thermal conductivity 
Prandtl number 
Heat flux 

Kin^atic viscosity 
Tank liquid height 
Tank liquid volume 
Gravitadcm acceleration 


TABLE 2 LIQUID OXYGEN TANK 


Svmbol 

Unit 

P 

psia 

To 

R 

D 

ft 

H 

ft 

V 

cuft 

d 

ft 

Q 

Cu ft/hr. 

P 

1/R 

P 

Ibm/cu ft 

Cp 

Btu/lbm. R. 

P 

lbm./ft hr. 

k 

Btu/fthr.R. 

Pr 


q” 

Btu/hr.sq ft 

V 

Sqft/hr. 

L 

ft 

VI 

Cuft 

g 

ft/sqhr 


Value Reference 

35 

140.5 

10 

23.83 

1589 

0.4 

1026.74 


0.02167 [ 2 ] 

74.99 [ 2 ] 

0.4 [ 2 ] 

0.687 [ 2 ] 

0.09685 [ 2 ] 

2.837 [ 2 ] 

29.681 

0.00916 [ 2 ] 

21.215 

1525.44 


417312.000 


L02 Tank 


Modified Grashof number 

Gr* = 6.691 X 10^^ 

Critical Reynolds number of L02 tank 
with selected D, H, and L 


Re = 4.4468 x 10^ 
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Inside Pipe 
Diameter 


(in.) 

Area (sq ft.) 

Velocity ft./s 

Reynolds No. 

Critical Re 

1 

5.45E-03 

5.23E+01 

1.71E+06 

4.45E+05 

2 

2.18E-02 

1.31E+01 

8.56E+05 

4.45E+05 

3 

4.91 E-02 

5.81E+00 

5.71E+05 

4.45E+05 

4 

8.73E-02 

3.27E+00 

4.28E+05 

4.45E+05 

5 

1.36E-01 

2.09E+00 

3.43E+05 

4.45E+05 

6 

1.96E-01 

1.45E+00 

2.85E+05 

4.45E+05 

7 

2.67E-01 

1.07E+00 

2.45E+05 

4.45E+05 

8 

3.49E-01 

8.17E-01 

2. 14E+05 

4.45E+05 

9 

4.42E-01 

6.46E-01 

1.90E+05 

4.45E+05 

10 

5.45E-01 

5.23E-01 

1.71E+05 

4.45E+05 

11 

6.60E-01 

4.32E-01 

1.56E+05 

4.45E+05 

12 

7.85E-01 

3.63E-01 

1.43E+05 

4.45E+05 
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